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Abstract 

The spin stiffness p s of the repulsive Hubbard model that occurs in the hy- 
drodynamic theory of antiferromagnetic spin waves is shown to be the same 
as the thermodynamically defined stiffness involved in twisting the order 
parameter. New expressions for p s are derived, which enable easier interpre- 
tation, and connections with superconducting weight and gauge invariance 
are discussed. 
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1 Introduction 



When a continuous symmetry is broken a Goldstone mode appears, which has 
to be considered as a slow variable and which necessitates an extension of the 
hydrodynamic equations Jj]]. The most famous example of this phenomenon is 
the case of superfluidity, where the Goldstone mode leads to second sound. The 
second sound velocity is governed by the superfTuid density p s , which is a measure 
of the stiffness of the superfluid order parameter against spatial variations. This 
situation has a perfect analogy in antiferromagnetic ordering, where the Goldstone 
mode is a spin wave and p s a spin stiffness. A third member in this family is 
superconductivity, where p s may be viewed as density of superconducting carriers. 
In the last example, one can relate p s also to the response to an electromagnetic 
field: it is inversely proportional to the square of the London penetration depth in 
the Meissner effect [||. 

In this paper, we analyze the antiferromagnetic ordering and the associated 
spin waves in the context of the Hubbard model. However, our considerations 
are applicable to a larger class of quantum lattice models; the Hubbard model 
merely serves as an example and an opportunity to give explicit formulae. The 
purpose of reconsidering this well-known theory is that it enables us to establish, 
on a microscopic basis, the general connection between spin-wave velocity and the 
stiffness p s [0. In doing so, we find new expressions for p s which reveal more clearly 
the nature of p s than does the rather formal definition in terms of a twist in the 
order parameter By using the Hubbard model as an example we can easily make 
the connection with the superconducting language, by means of a transformation 
between repulsive and attractive Hubbard models. This gives us the opportunity 
to clarify some misunderstanding in the literature concerning gauge invariance. In 
general, the spirit of the paper is an articulation for the Hubbard model of the 
analysis by Forster |]J of the connection between broken symmetry, correlation 
functions and hydrodynamics. 

We have organised the paper as follows: we start with a brief discussion of the 
tools needed for the analysis. Then we show how p s enters in the (hydrodynamic) 
equations for the spin waves. In section 4, we work out the (thermodynamic) 
definition of p s in terms of a twist in the order parameter and subsequently we 
demonstrate in section 5 that the hydrodynamic p s is indeed the spin stiffness. In 
section 6, we draw the analogy with p s as the response function between the current 
and the inducing gauge field as it appears in the theory of superconductivity. In 
section 7, we comment on formal manipulations with the expression for p s and derive 
new expressions which are more transparant. We close with a brief discussion. 

2 Preliminaries 
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2.1 The Hubbard model 



The Hubbard model is represented by the hamiltonian ||: 

3,S,a j 

where Cj} are electron annihilation (creation) operators for site j with spin a 
(a =t, I corresponds to +1,-1, respectively, if it does not occur as an index). 
Neighboring sites of site j are denoted by j + 5. We consider a bipartite lattice, 
such that hops from site j to site j + 5 are always from one sublattice to the other, t 
is the exchange- or hopping integral, U is the on-site interaction between electrons 
of opposite spin, and rij a = Cj a Cj a is the occupation number operator. 

Tih has a rich set of symmetries. For our purpose it suffices to mention the 
magnetic symmetry: 

[?4,M]=0, (2.2) 
where M is the total magnetization operator: 

M = $>, (2.3) 

3 

with Sj the local magnetization operator: 

Sj = lE C^CT^'CjV , (2.4) 

with a the set of the three Pauli matrices. 



2.2 Averages and operator inner products 

Averages are based on a (grand) canonical ensemble described by a hamiltonian TC. 
The partition function gives the free energy F according to: 

F = — i lntr e _/3H , (2.5) 
P 

where (3 = 1/ksT. Averages are taken with respect to the canonical weight: 

(A) =e pF tie- pn A. (2.6) 
We will need an inner product in operator space, for which we take || |l[: 

(A, 13) = ^J Q P d\{(A^- XH Be m ) - (A*)(B)} . (2.7) 

In addition to having all the properties of an inner product this definition has a few 
more desirable aspects. If the hamiltonian H is perturbed by adding the operator 
STi, the linear response of an operator J is given by: 

5(J) = (J) H+ sn ~ (J) = -P (SH\ j) . (2.8) 
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The change in the free energy is to second order in 5H.: 



5F = (5H) - - (5H\ 5H) . (2.9) 
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The last two equations are derived in the Appendix. 

Another property of the inner product, that we will use, concerns the Heisenberg 
operators: 

A(t) =e itn / h Ae- un/h (2.10) 

and their time derivatives 

A(t) = % - [H, A(t)] . (2.11) 

Using the invariance of the trace under cyclic permutation of the operators one 
easily proves the relation: 

(B,A(t)) = j K ([&,A(t)]). (2.12) 
2.3 Hydrodynamic equations 

Linear hydrodynamic equations are obtained by first identifying the set of hydro- 
dynamic variables Aj JI[ . They span a hydrodynamic subspace (using the operator 
inner product defined above). Projection on the hydrodynamic subspace is achieved 
by the projection operators: 

PO = '£A j P jm (,A m ,0) , (2.13) 

j,m 

where the matrix P is the inverse of the matrix of susceptibilities (Ai,Aj): 

Pjm {A m ,A k ) = 8 jk , (2.14) 

m 

where all indices run through the hydrodynamic subspace. 

The dissipation-free hydrodynamic equations result from the projections on the 
hydrodynamic subspace of the time derivatives: 

P^ = XMi n *» (2-15) 
j 

with Qj£ given by 

n j t = y EP jm (A m ,Ai) (2.i6) 

m 

The hydrodynamic equations then read: 

^W*)>=E<^(*)>%- (2-17) 
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3 Spin stiffness and hydrodynamic spin waves 



The slow or hydrodynamic modes derive from the conserved quantities and the 
Goldstone modes of a broken symmetry [|TJ. Since the total magnetization is con- 
served, the fourier components 

M(k) = eik ' rJs i (3- 1 ) 

i 

are slow modes for small k. The symmetry breaking that we consider is antiferro- 
magnetic ordering, which occurs in the Hubbard model at low temperatures near 
half-filling. The staggered magnetization 

N = ]T(-iy Sj (3.2) 

j 

then acquires a non-zero average. (— l)- 7 equals +1 on one sublattice and —1 on 
the other sublattice. N can be seen as a fourier component M(Q), where the wave 
vector Q has the property: 

e iQ-r, = (33) 

The fourier components N(k) = M(k+Q) are the slow modes for small k associated 
with the antiferromagnetic symmetry breaking. We assume that the system orders 
antiferromagnetically in the x-direction, i.e. 

(N x ) = M st ^ . (3.4) 

The spin waves arise from an interplay of the small-k components of M(k) and 
N(k). Which components play a role can be seen from the general commutation 
relation 

M a (k), M /3 (k / )] = i M 7 (k + k') , (3.5) 

with a, /3,7 a cyclic permutation of x,y,z. Taking 7 = x, we see that a non-zero 
average of N x couples M z (k) to N y (k) and M y (k) to N z (k). The two pairs are 
equivalent and we focus on the first couple. In principle, one should write down the 
hydrodynamic equations in the combined 6-dimensional space of M(k) and N(k), 
but symmetry considerations permit to restrict the equations to the 2x2 subspace 
of M z (k) and N y (k). So we obtain equations of the form |l|, [/J: 

<)(1[ ' lk - t ">l = (M z (k,t))Q zz + (N y (k,t))Q yz , (3.6) 



m 

d(N y (k,t)) 
dt 



(M z (k,t))n zy + (N y (k,t))n yy (3.7) 



The VLj£ are calculated with (|2.16|) . The entry (A m} -A.e) is simple, as it can be 



calculated with (|2.12j ) and (|3.5|) : 

(M z (k),M z (k)) = (M z (k),N y (k))=M st /(3h (3.8) 

(N y (k),M z (k)) = -M^/ffi (N y (k),Ny(k)) =0 . (3.9) 
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The matrix elements Pj m are obtained from the inverse of (A m ,Ak), which in turn 
is interpreted as the response to a perturbation as in fl2.8|) . So we consider first a 
perturbation: 

5H = -£M z (k), (3.10) 
with a small amplitude £. It will induce: 

5(M z (k)) =/?£(M 2 (k),M 2 (k)) =iV X± (k)e , (3.11) 

where iV is the number of sites. In the second equality we introduce a susceptibility 
X±, where the subscript is a reminder that it concerns a susceptibility perpendicular 
to the order parameter direction. We may drop the k-dependence since we are 
interested in small k and x_i_(k) is expected to be regular for small k. 

For symmetry reasons the perturbation ( ft .10] ) will give no effect 5(N y (k)). So, 

(M 2 (k), N y (k)) = (N y (k), M z (k))* = . (3.12) 

The perturbation, 

6H = -£N v Ql) , (3.13) 

requires a different treatment. It tends to turn over the spontaneous order from 
the x-direction to the y-direction. For small k the response will diverge and we will 
determine the singular behavior in section 5. For the moment, we use the result, 

M 2 

(N y (k), N y (k)) = j^j- 2 (1 + 0{k)) k - . (3.14) 



One can see ( 3.14 ) as the hydrodynamic definition of the spin stiffness p s . 
Now all the ingredients for Qje are present and we find: 

tt zz = tt zv = — — (3.15) 

n yz = --jjfY n vy = ° ( 3 - 16 ) 



The hydrodynamic equations (ft.6| ) -( pl7|) have a solution: 

(M z (k,t)) = (M z (k,0))e ±tckt (N y (k,t)) = (N y (k,0))e ±icht , (3.17) 
with the spin wave velocity c given by: 

c=Lp±, (3.18) 
n V X± 

a relation due to Halperin and Hohenberg || . 
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4 Spin stiffness and twisted order parameter 



The more fundamental definition of the spin stiffness p s relates it to the increase 
in the free energy due to imposing a twist in the boundary condition on the order 
parameter |4[]. E.g., one can have the order parameter point in the ^-direction at 
one end of the system and in the ^/-direction at the other end. For a continuous 
symmetry the twist in the order parameter will spread itself equally over the system. 
So we imagine that the order parameter slowly rotates around the z-axis: 

(sj) =m st cos ((q+Q) -Tj) (sj) = TO* sin ((q + Q) • r,-) , (4.1) 

where q is the pitch of the twist and m st = M st /N. The two equations ( |4.1| ) can 
be combined into: 

= m st exp (i (q + Q) ■ r d ) , (4.2) 

with 

4 = s] + is] = c^c n , (4.3) 

the spin-raising operator. The condition ( [4. 2D leads to an increment in the free 
energy which for small q can be written as: 

F(q) = F(0) + \Np s q 2 + ■■■ . (4.4) 

The term linear in q is absent due to inversion symmetry. Formula ( (4.4|) has to be 
considered as the thermodynamic definition of the spin stiffness p s against a twist. 

This formal definition does not lead to a simple calculational scheme since it 
involves the computation of the free energy of a spatially inhomogeneous system. 
We can however map the system with the condition ( |4.2j ) on a different system with 
an easier constraint by an unitary transformation U of the form: 

U = e iq ^ T *"i . (4.5) 

U transforms H, and the density matrix p as: 

W = WHU ] p' = UpU^, (4.6) 

such that the partition function is invariant: 

e -^ = t re-^ = tre-^' . (4.7) 

Averages in the transformed system (based on p') are denoted as (•■■)' and are 
related to the original averages (• • •) by 

(A)' = tr p'A = tr UpU ] A = trpU^AU = {A'} , (4.8) 

with 

A' = WAU . (4.9) 

U locally rotates the spins around the z-axis over an angle q-r^. Using the properties 
of spin operators, we have: 

s f = uh+u = c ;,, r .s; , (4.io) 
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such that (|4.2|) changes into: 

(sp'=(sf)=m st e^ , (4.11) 

which puts the staggered magnetization in the new system everywhere in the re- 
direction. Thus the constraint is simple and homogeneous at the expense of chang- 
ing the hamiltonian from H to H'. It is however easy to compute the new Hubbard 
hamiltonian H' h by making use of the properties of the Cj a under the transformation 
U: 

Wc^W^e^V. (4.12) 

Using 0, H' h is: 

K = ~t E ^ rs/2 c] + s,^ + Uj2n n n jb (4.13) 

j,8,<r 3 

where = rj + g — Xj is a vector connecting nearest neighbors. Note that the new 
hamiltonian is again spatially homogeneous and of the Hubbard form with a com- 
plex hopping integral that depends on the pitch of the twist. So we must evaluate 
the free energy of H' h under the condition that the order parameter points in the 
^-direction. In, e.g. the mean-field approximation one can directly determine F(q) 
and then expand F(q) to obtain p s ||. Concrete calculations of p s starting from 
the thermodynamic definition (using series expansions) can be found in Refs.0 and 
[p~0[1 for the 2D Heisenberg antiferromagnet and Hubbard model, respectively. 
In general we may evaluate F(q) by expanding H' h in powers of iq: 

H' h = H h + iq ■ Hi + izq • % ■ iq + • • • , (4.14) 

with 

#i = ~E r 5 ac] +Sta c ja = NJl (4.15) 

n 2 = -j E r s r sc] +5 ,a c j« ■ ( 4 - 16 ) 

3,$& 

We note that Hi is the current of the ^-component of the spin and that 7i 2 is (for 
cubic symmetry) a diagonal tensor with the kinetic energy on the diagonal. As q 
is small, the free energy increment due to q can be calculated using perturbation 
theory using ( p.9| ): 

F(q) - F(0) = |iq • { (%) - (3 (h\, Hi) } • *q (4.17) 

Here we use that (Hi) = (no current in the unperturbed state) and we note that 
H[ = —Hi. So we obtain for p s the expression: 



\ {(H 2 )* + P (HiMiY) ■ (4.18) 



Ps= N 

For cubic symmetry p s reduces to p s = p s I. We stress, by using superscripts x, that 
in ( |4.18| ) the averages are with respect to an ensemble with the order parameter 



pointing in the x-direction (and with Hh as hamiltonian). 
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5 Equivalence of the two definitions for spin stiff- 
ness 

We now prove that the hydro dynamic definition Q3.14 ) and the thermodynamic 
definition ( |4.4| ) are equivalent. We do this by considering the response to the 
perturbation: 

«H = -E(-lft<, (5.1) 

i 

where the £j are small and smoothly varying in space, e.g. of the form: 

& = £ cos (q ■ Tj) . (5.2) 

The perturbation is added to the hamiltonian: 

H = H h - CN X , (5.3) 

where we have introduced a small symmetry-breaking field £ X N X to guarantee that 
the reference system has its order parameter pointing in the x-direction. 

Because of the intrinsic stiffness of the system the response to the small per- 
turbation ( |5.1|) will be a local rotation of the spins around the z-axis over a small 
angle iff 

(s*) = (-1) 3 m st cos (pj , (s y j) = (-1) 3 m st sin cpj . (5.4) 

where <pj is spatially smooth. The spatial variation (|5.2j) will induce a similar wave 
pattern for the (pf 

(fj = A cos (q ■ Tj) , (5.5) 

with a small amplitude A proportional to £. Determining this proportionality is 
the main goal of this section, as it leads directly to the singular behavior of the 
inner product ( |3.14| ). To see this we write 57i with (5T2) for the £j as: 

5n = -^[Ny( q ) + Ny(- q )] . (5.6) 

The anticipated response ( |5.4| ) gives for the ^-component of N(q): 

h(N y (q) + N y (-q)) = m st £ cos(q • Tj ) sin^ . (5.7) 

3 

Inserting ( |5.5|) and using that the angles <pj are small yields: 

(iV%) + jV»(-q)) = 2m st A^ cos 2 (q • Tj) = M st A . (5.8) 

3 

On the other hand, linear response theory (see ( |2.8| )) gives: 

(W(q) + W(-q)> = ^ (W(q) + W(-q), N*(q) + ^(-q)) . (5.9) 



The modes q and — q have no overlap, so combining ( |5.8j ) and ( |5.9p leads to: 

M st A = /3£(iV%),iV%)) . (5.10) 
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Thus the computation of A gives the inner product which must lead to the hy- 
drodynamic definition of p s , ( |3.14 ). We calculate the ipj by first determining the 
increase in the free energy for a set of arbitrary ipj and then optimize this free 
energy with respect to the tpj] this yields the cpj resulting from the £j. As in the 
previous section, we gauge away the ipj by an unitary transformation: 



The averages in the new system are: 

«)' = (Z/sJW) = (-l)%i st 



(Uh y M) =0 



(5.11) 

(5.12) 
(5.13) 



such that in the new system the order parameter points in the x-direction every- 
where. This new constraint has to be combined with the transformed hamiltonian: 



H'=U\H h - CN X - J2 u] ■ 



(5.14) 



We evaluate the free energy by perturbation theory. Therefore we decompose "H' 
into: 

W = H h + H a + H h , (5.15) 

where Tih is taken as the unperturbed hamiltonian together with the constraints 
( 5.12|) - (|5.13|) on the order parameter. As in the previous section, such averages are 
denoted by (• • -) x . The other two parts in ( |5.15| ) are given by: 

= un h u^-n h 



-u 



(5.16) 
(5.17) 



These parts will be treated as perturbations on Ti h . can be expanded for small 
angles as: 



Y,'Hjrj -X <PiHij<pj + -- 



h3 



where it follows formally from ( 5.11| ) that 



s*,H h 



For the Hubbard hamiltonian (cf. Q4.13) ) one can explicitly transform Tih. 



urt h w 



-*£« 



(5.18) 

(5.19) 
(5.20) 

(5.21) 



which can then be expanded as in ( |5.18|) , yielding explicit expressions for Tij and 
Tiij. The free energy associated with 7i a can be written as: 



(5.22) 



hi 
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with 

K ij = -{n ij ) x -i3{n i ,n j ) x . (5.23) 

This part of the free energy increase can be directly related to p s as defined in 
the previous section. Kij is translationally invariant, i.e., a function of r« — r^. 
Moreover, we have: 

K ij = K ji , £i^ = 0, ( 5 ' 24 ) 

3 

as a constant rotation angle (pj will not lead to an increase in the free energy. 
Angles with an uniform gradient, 

(Pj = q- Tj , (5.25) 

have been discussed in the previous section and lead to p s . So: 

£ (q ■ Ti) (q • r,-) = N q_ ■ % ■ q . (5.26) 

i,3 

We will need the fourier transform of K^, which can be evaluated with Q5.24 ) and 
(|5T26| ) as: 

J] ^ e^'^ = N q ■ % ■ q q - . (5.27) 

The part TCh is found as: 
H h = -J2 { tr cos ^ + sin ^] sj + [-£ a sin ^ + & cos ^] sj} , (5.28) 

3 

and yields: 

SF h = (H b ) x = -m st £ (£* cos ^ + sin ^) . (5.29) 

3 

Thus the total increase in free energy is: 

SF = 5F a + 5F h = i <-Pi K ij<-Pj ~ m st 0T cos <pj + ^ sin ipj) . (5.30) 

i,3 3 

Now the follow from the stationarity of 5F, which leads (for small <pj) to the 
linear equations: 

Y + m st £ x (pj = €j m st . (5.31) 

Due to the translational invariance of one finds that indeed ( |5.2| ) and ( |5.5| ) form 
a solution with A and £ related by: 

(K(q)+m a £ B )A = m 8k Z , (5.32) 

where K(q) follows from ( |5.27| ): 

K(q)=q-%-q q^O. (5.33) 

Now we can safely let £ x — > and one sees that ( |5. 10| ) , ( [5.32 ) and (|5.33|) lead to: 

(iV»(q),JV»(q)) = — ^ q-0, (5.34) 

pNq • p s - q 

which is a slight generalisation of ( |3.14f ) to systems with a lower symmetry than 
the cubic symmetry. Since indeed the same ~p s appears as in the previous section, 
we have demonstrated the equivalence of the thermodynamic and hydrodynamic 
definitions of the spin stiffness. 
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6 Response to a Gauge Field 



A third aspect of p s , which is unusual in the context of antiferromagnetic order, is 
its role in the response to a gauge field. This role is familiar in superconductivity, 
where the presence of p s leads to special electromagnetic behavior in a symmetry- 
broken state: the Meissner effect. In a lattice system the electromagnetic field is 
introduced by the so-called Peierls substitution [[□]]. It amounts to the replacement, 



Cj a Cj a exp 



-le 



P dr- A(r)/ft 
Jo 



(6.1) 



where e is the electric charge and A(r) the vector potential. Executing this re- 
placement in the Hubbard hamiltonian changes it into: 

We have introduced the abbreviation: 

[ rj+6 dr • A(r) = A j+5/2 ■ r 5 , (6.3) 

such that in the lattice version the vector potential Aj +< 5/ 2 is associated with the 
link between the sites j and j + 5. By introducing the vector potential in this 
way one has to change the expression for the electric current concomittantly. The 
current Jj+s/2 is also associated with the link \j,j + S] and it should obey the lattice 
version of charge conservation, 

^ + E W ■ * = , (6-4) 
with the charge density given by: 

Pi = e E n i° ■ ( 6 - 5 ) 
Using the equation of motion for pj one finds: 

W = ^ £ {c] +5>a c ja e ieA ^rr s /n _ c ) a c 3+ ^-^ + ^) . ( 6 .6) 

When the system becomes superconducting the electric current J becomes propor- 
tional to A (the London equation): 



rn 



(6.7) 



where n s is the density of the superconducting carriers and m is their mass. Super- 
conductivity occurs in the negative-f/ Hubbard model and it is the counterpart of 
the antiferromagnetism in the positive-f/ Hubbard model that we have been con- 
sidering so far. The two can be related to each other by the transformation (see 
e.g. 0): 

C ; T = c iT , c' n = (-iy c] l . (6.8) 
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The transformation (|6.8| ) transforms the charge density, 

p'j = e H n 'j° = e E an i° = 2e s j ' ( 6 - 9 ) 



into the z-component of the magnetization. The hamiltonian 1~La from (|6.2| ) is 
transformed into: 

h 'a = ~t E c) + ^c^ A ^/^l h + UY, • (6-10) 

j,S,a 3 

Note the appearance of a factor a in the phase factor. Comparing this expression 
with ( |5.21| ), we see that they are equivalent for: 

A i+5/2 ■ r<5 = (ip j+s - tpj) y & ■ (6.H) 

The charge conservation law transforms into a conservation law for the z-component 
of the magnetization: 

Qs z . 

+E J^/2-r 5 = 0, (6.12) 

with the current Jj+s/2 given by: 

The equivalent of the London equation is the response of J z to A. We have two 
contributions: one because J z depends on A and one because the hamiltonian 
(|6.10| ) depends on A. Expansion of J z gives: 



Jj+s/2 = J jW + T" J1 i+^/2 • A j+5/2 H , (6.14) 



2e = 



'j+8/2 - u j+8/2 -r ^ 

with 



2h 

^3+6/2 = -^-r 5 r 5 J2 { c }+8,aCja + c^cj+s^) ■ (6.16) 

Secondly, we have a contribution which results from the expansion ( 6.1(J| ) for small 
A: 



iet 

3,8,<r 



PHa = E 4+8,a C ^ r 8^ ■ A j+S/2 ■ (6.17) 



We note that Hi and 7^ 2 as found in ( f4.15|) - ([4.16|) are related to J° and J 1 (Note 
that the two terms in the sum over cr in ( |6.15|) and Q3.16j ) contribute equally): 

n, = (6-i8) 

% 3,8 

^2 = 3+8/2 (6.19) 

3,8 
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The London equation (|6.7|) only holds for sufficiently slowly varying Aj + s/ 2 , i-e. 
the spatial variations must be small over a coherence length. So we may take A 
constant and we then also have: 



5H A 



2ie 



Hi- A 



The total magnetization current is defined as: 

1 



3,S 



2N ^ 1+5/2 



From linear response theory we obtain for homogeneous A (cf. 
(2 z )h+8h 



1 EtM'A-^ ^,e j ; +( v- 

3,6 11 \ j,S 

2e r ,= 



2iV 



{(H 2 )+P(Hi,Hi)}-A. 



Comparing this with the expression (|4.18 ) for p s , we may write: 

(J Z )h+5H : 



2e 

n 



T2 Ps ■ A ■ 



(6.20) 



(6-21) 



(6.22) 
(6.23) 

(6.24) 



Therefore indeed the same p s appears as before and we see that the superconduct- 
ing weight in the response to a gauge field is also equivalent to the spin stiffness 
associated with a twist in the order parameter. Very recently, a discussion of the 



effect of a spin-dependent gauge field, as used in this section, was given in Ref. |L2 
in the framework of Fermi liquid theory. 



7 Alternative forms for p s 

In this section, we discuss some formal manipulations with the expression for p B . 
These are included because they involve some subtleties due to the symmetry break- 
ing. The manipulations lead to expressions which are more general or easier to 
interpret. The formula (|4.18|) for p s together with the definitions (|4.15|) and (|4.16| ) 
refers explicitly to the Hubbard model. We can obtain a more general expression 
by expanding U, as given by (fL5[) in powers of iq: 



U 



1 + iq • d z + 



(7.1) 



with d z given by: 



r = E 



(7.2) 



Inserting ( [7.1| ) in the definition of 7i' yields a similar expansion as ( |4.14j ) with: 



n 1= r,n h 



(7.3) 
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Using this in the expression (|4.18| ) gives for p s : 



Ps 



(7.4) 



This expression makes no reference to the Hubbard model and is therefore gener- 
ally valid for lattice hamiltonians. We have augmented the average and the inner 
product with a superscript x to reflect that the averages are taken in a system 
where the order parameter points in the x-direction. Such a warning signal is not 
superfluous as the following "derivation" may show. Suppose we use the relation, 



■&* = - 
a 

to write for the second term in (|7.4[) : 



h 

% 



and use the relation (|2.12|) to transform it into: 



r) z 



$ Z ,H ) . 



(7.5) 



(7.6) 



(7.7) 



Then we would discover that the two contributions in ( |7.4|) exactly compensate 
each other! The error in this reasoning is that we have ignored the fact that the 
averages have to be taken in a symmetry-broken state. However, relation Q2.12| ) is 
based on a cyclic rotation of the operators in a trace, which is only permitted when 
the trace is taken over the whole Hilbert space. 

We can elucidate this point further by implementing the constraint on the order 
parameter by including a symmetry-breaking term in the hamiltonian as in (|5.3|) . 
Then we would have for (|7.5|): 



■& Z ,H +C 



The extra term looks innocent, because £ x is vanishingly small, but it is not, as we 
shall see. Using (|7.8|) , gives for the second term in ( [7.4|) : 



([$ z ,H h ],[& z ,H h 



if 



(7.9) 

Using ( |7.8| ) again in the first entry of the second term of ( |7.9| ), we obtain two 
contributions, one proportional to and one to {l; x ) 2 - The former reads with the 
use of 



e 



(3h z 



J2 r i r i ! "'s' • 



(7.10) 
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The lattice sum can take arbitrarily large values and thus one cannot rely on the 
smallness of £ 2 to ignore its contribution. The term proportional to (£, x ) 2 is equally 
difficult to interpret. 

We may however use these types of operation to give the expression for p s yet 
another form, and thereby bring this subtle point into focus. First, we relate the 
second part of (|7.4J) to a current-current inner product: 



([r,H h ],[# z ,H h 



h 2 N 2 (Jt, J 



hi °h) 



with J I defined in ( |4.15|) . Next, we use the identity, 



r) z 



d 1 



■& x ,H h 



(7.11) 



(7.12) 



which is a direct consequence of the magnetic isotropy and which is easily proven for 
the Hubbard model by direct evaluation. Now for the right-hand side of ( |7.12| ) we 
can apply the train of arguments (|7.5|) - ([7.7|) everywhere replacing z by x, since d x 
commutes with the symmetry-breaking term. So (|7.7| ) holds for the x-component. 
This permits us to write the expression for p s in the form: 



Ps 



^ 2 iV{(j|,jl)-(j^,J^)} 



(7.13) 



where we have used the coordinate-free notation: parallel and perpendicular are 
to be understood with respect to the orientation of the order parameter. This 
expression most clearly shows that p s is induced by symmetry breaking: without 
symmetry breaking the distinction between parallel and perpendicular disappears 
and p s vanishes. In a mean-field (or: BCS) approximation one finds at T = only 
the parallel contribution. For increasing T a gradual compensation occurs between 
the two terms, which becomes complete at T = T c . 

Since ( [7.12j) relies on the magnetic isotropy of the model, the expression (|7.13 ) 
is no longer valid when the isotropy is broken by an external magnetic field in the 
^-direction. Then ( |7.4j ) still holds and also the Halperin-Hohenberg relation ( |3. 18| ) 
remains valid. Of course both p s and x± are affected by the presence of such a 
magnetic field, as can e.g. be seen from a mean-field treatment of these quantities 
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As a final comment on this genre of expressions we discuss the Bogoliubov 
inequality which is sometimes used to make ( |3.14| ) plausible M. The general form 
reads: 

(A,A)(B,B) >| (A,B) | 2 . (7.14) 

Take A = N y (k) and B = M 2 (k) such that (A, A) is the desired inner product. 
The right-hand side of (|7.14j) is then given by (|3.9|) . One might think that (B,B) 
is associated with the inner product (J 2 , J 2 ), since for small k one has: 



M 2 (k) ~ - • ti z ,H 
n 1 



ziVk-J 2 . (7.15) 
This is not correct as a more precise analysis shows. First we use ( [2.12| ) to write: 



(M 2 (k),M 2 (k 



i 

]5h { 



M 2 (-k),M 2 (k) ) . 



(7.16) 
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Then we expand the first entry as in (7.15 ), but taking the symmetry-breaking 
term explicitly into account: 



M*(-k) 



V(A^(-k) - N y ) . 
n 



(7.17) 



The small symmetry-breaking term gives: 







■([N y (-k)-N v ,M z (k)\) 



£ X M S 



St. 



Then we obtain for small k for (|7.16|) : 
(M*(k),Af*(k)) = ^j{-k. 



>-k + fM st } 



Substituting all of this in ( [7. 14j ) we have: 
{N y {k),N y {k)) > 



Ml 



/3(-k-([^, [&,H h ]})-k + ^M s 



(7.18) 



(7.19) 



(7.20) 



This result should be compared to ( |5.34|) ; substituting for p s the expression (|4.18|) , 
with ( |7.3|) for TCi and H2, one sees that the Bogoliubov inequality above only 
involves the double commutator, or, equivalently, the inner product of parallel 
currents, and not the inner product of perpendicular currents. Note that (|7.15f) 
seems to suggest that the inner product of perpendicular currents would be involved. 
This is not the case and generally we are indeed faced with an inequality (the 
symmetry-breaking term is innocent in this respect). 



8 Discussion 

In the foregoing sections we have discussed the role of p s : 

• as a parameter in the hydrodynamic (spin wave) equations, 

• as the parameter determining the increase in the free energy due to long- 
wavelength variations in the phase of the order parameter, 

• as the proportionality parameter between a perturbing gauge field and the 
induced current. 

By explicit calculation we have shown that it is the same p s entering in all these 
aspects. We have placed the discussion in the context of the Hubbard model which 
allows more explicit expressions. As we have indicated, the expressions can easily 
be extended to more general lattice models. 

In the case of magnetic isotropy we can write p s as the difference of two current- 
current inner products containing the spin currents parallel and perpendicular to 
the order parameter. The Bogoliubov inequality only involves the parallel current. 
In the BCS approximation it becomes an equality only at T = 0. 
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A final comment we want to make concerns the issue of gauge invar iance. The 
London equation flb.7| ) can be extended to the frequency-wavevector domain as: 

(J(k,w))=^(k,w).A(k,w), (8.1) 

where A(k, ui) is the double fourier transform of A(r, t) with respect to space and 
time and x(k, oS) a generalized susceptibility (which has the role p s had before). The 
limit of slow variations in space and time is delicate ||. In the limit of first k — > 
and then u — > one has electric response and x leads to the Drude weight. In the 
opposite order, first u — * and then k — > 0, one has magnetic response leading to a 
no n- vanishing superconducting weight, which is considered here. Electromagnetic 
gauge invariance requires that a purely longitudinal A(k, uj) = ka(k, uj) has no 
response, since it can be gauged away. Thus: 

f(k,^)-k = 0. (8.2) 

We have been considering static response throughout, so we have taken the limit 
uj — > first and we have not been worrying about the k-dependence of the spin 



stiffness p s . Thus the expressions given for p s do not obey the condition (p72|). 
The origin of this paradox is not to be traced to shortcomings of e.g. the BCS- 
approximation, as is sometimes done in the idea that vertex corrections will restore 
gauge invariance 0. The answer comes from the change in gauge invariance in 
the symmetry-broken state. In section 6, we saw that gauge fields are similar to 
imposed twists on the phase of the order parameter. The theory is required to be 
invariant under the combined gauge transformation: 

A 'j+s/2 ■ r s = A i+5/2 ■ r s + Xj+s - Xj (8-3) 
2e 



i = h + irxj, (8-4) 



where <f)j is the phase of the order parameter and Xj an arbitrary function. Corre- 



spondingly the London equation (|6.24fj should be formulated in a gauge-invariant 
way as: 

(JW = -p7V U- . (8.5) 

Then the requirement of gauge invariance is automatically fulfilled and it imposes 
no further requirements on p s . 
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Appendix 



In this appendix we derive expressions for the change in the average of an operator 
J and the change in free energy upon perturbing the hamiltonian 7i by adding an 
operator 57i. These expressions, up to first and second order in 5H, respectively, 
are given in equations ( |2.8| ) and ( |2.9| ) in terms of the inner product defined by ( |2.7| ) 
and are repeated here: 

S(J) = (J) n+ 8H-(J) = -P(5H\j), (A.l) 
l ntre -m+SH) _ lntre -/wl = (SH) - - (5H\5H) . (A.2) 



Unless explicitly indicated averages are taken with respect to the operator e" 
(see©). 

The derivation of both equations proceeds by way of the operator identity: 



e -(3(H+SH) = e -tm 



f d\e~^ n 6ne-^ n+5H \ (A.3) 



which is proven by noting that both left- and right-hand side are equal for (3 = 
and have identical derivatives with respect to (3. 

Equation ( |A.1| ) is now derived by inserting ( |A.3| ) into: 

{J)h+sh - tre _ m+m) , (A.4J 

expanding to first order in 57i and using the following result (which is easily proved 
by cyclic permutation in the trace): 

[ P d\(5H(\)) = f d\(e xn 5He- m ) = (3(8H) . (A.5) 
Jo Jo 

To derive ( |A.2j ), we start by iterating ( A.3|) once and expanding to second order 
in 5H: 

tTe -m+m = tr e -w L _ d x(5n(X)) + £ dA £ d\ r (5H(\)5H(\')) + ■ ■ ^ . 

(A.6) 

Before taking the logarithm in order to compute 5F it is useful to prove the identity: 

[ P d\ t d\'(6H(X)6H(\')) = ^ f d\(6H(\)5H) . (A.7) 
Jo Jo 2 Jo 

The proof proceeds by a sequence of changes of variable in the integral, introducing 
the notation: 

f(r) = (5H(t)5H) , (A.8) 

and using: 

/C9-r) = /(r) . (A.9) 
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Equation ( |A.9|) is again proved by cyclical permutations in the trace. Explicitly: 



/ dA [ X d\'(5H(\)5H(\')) = / dA [ X d\'(5H(\ — \')5H) 
Jo Jo ' Jo Jo 

/ dA / dr/(r)=/ dr/ dA/(r) = / dr (/3 - r)/(r) = 

JO JO JO Jr JO 

1 f/3 



dr r f(P - r) = i f dr (/? - r + r)/(r) = 2 T dr (5H(t)5H) 
z Jo 2 JO 



Now equation ( |A.2j ) is easily derived by inserting ( |A.7|) and ( |A.5|) into ( |A.6| ), taking 
the logarithm, expanding to second order in STi, and recognizing in the final result 
the inner product (|2.7|). 
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